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ABSTRACT 

Using the data obtained with the Spitzer Space telescope as part of the Surveying the Agents 
of a Galaxy's Evolution (SAGE) legacy survey, we have studied the variations of the dust com- 
position and abundance across the Large Magellanic Cloud (LMC). Such variations are ex- 
pected, as the explosive events which have lead to the formation of the many HI shells observed 
should have affected the dust properties. Using a model and comparing with a reference spectral 
energy distribution from our Galaxy, we deduce the relative abundance variations of small dust 
grains across the LMC. We examined the infrared color ratios as well as the relative abundances 
of very small grains (VSGs) and polycyclic aromatic hydrocarbons (PAHs) relative to the big 
grain (BG) abundance. Results show that each dust component could have different origins or 
evolution in the interstellar medium (ISM). The VSG abundance traces the star formation ac- 
tivity and could result from shattering of larger grains, whereas the PAH abundance increases 
around molecular clouds as well as in the stellar bar, where they could have been injected into 
the ISM during mass loss from old stars. 
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1. Introduction 



The Large Magellanic Cloud (LMC) is the nearest galaxy external to the Milky Way (MW), located 
at a distance of about 55 kp c, and is oriented with an advantageous almost face-on viewing angle (35°, 
van der Marel & Cioni 11200 ih . The LMC is a dwarf galaxy orbiting around the MW, containing more than 
30 billion stars. As opposed to the situation in the Galaxy, where studies of the infrared emission from dust 
in the interstellar medium (ISM) suffers from significant confusion along the line of sight, the favorable 
viewing angle for the LMC offers a direct view to the processes taking place in the diff use ISM. The LMC 
has a lower metallicity tha n that of our Galaxy, with estimated values close to 0.3-0.5 Z© (|Westerlund 1 ll997h 



or 0.25 Zq (iDufour II 19841) ; specifically for carbon and oxygen, Z[^j^^ 



= 0.28 Z^ and Z^^^ = 0.54 Zg. 



Extinction curv es in several regions of the LMC (IClayton & Martin II 19851 : ISauvage & Vigroux II 1991 



Gordon et al. Il2003h also differ than those measured in the Galaxy. For example, the ultraviolet (UV) bump 
at 2175 A is weaker in the LMC than in the MW. Moreover, the far-ultraviolet (FUV) rise is steeper in 
the LMC, indicating that grains responsible for this feature could be more abundant. However, the LMC 
and M W extinction curves are similar in the visible (VIS) and the near-infrared (NIR). The Nhi/Av ratio 
values (IGordon et al. II2003I : ICox et al. Il2006r) ind icate a gas/dust ratio ranging between two and four times 
the Galactic value, depending on the line of sight. Isakon et alT (2006) analyzed the LMC far-infrared (FIR) 
emission using DIRBE data and showed the existence of some pixels with very bright emission at 240 yum 
which could possibly be interpreted as cold dust. Emission in DIRBE bands at 12 and 25 /im normalized 
to FIR brightness was found to be relatively lower in the LMC than in our Galaxy. These results can be 
interpreted as a polycyclic aromatic hydrocarbon (PAH) and very small grain (VSG) deficit in the LMC. 
Sauvage & Vigroux I (Il99lh studied the Infrared Astronomical Satellite (IRAS) 1 1 2/1 2.5 versu s I(,n/Ii nn ratio 



in the LMC and SMC and compared them with observations of several galaxies by iHelou I (Il986h . They 
showed that there is a decrease of the llfim emission relative to the 60 and 100 /im emission associated 
with the decrease of the metallicity in the LMC and SMC compared to the MW. They proposed the low 
metallicity in the LMC and SMC could limit small grain formation, concluding that the metallicity could 
act preferentially on the abundance of one grain component. 



Bernard et al. 



(l2008h pe rformed a global preliminary analysis of the entire LMC using new Spitzer 



images (Meix ner et al. 1120061) covering the entire LMC with much higher angular resolution than the previ- 
ous IRAS or DIRBE data. They produced a temperature map for the overall LMC and found temperature 
variations in the range 12.1 - 34.7 K. They studied three individual regions, determining and modeling their 
spectral energy distribution (SEDs). They showed that the relative abundance of the dust components is 
close to that in our Galaxy, in contrast to the previous IRAS and DIRBE studies. They also found that the 
PAH and VSG abundances are higher in the molecular phase than in the neutral phase for a few selected 
regions. They explained this result as an increase of the two types of grains in unresolved halos surrounding 
molecular clouds. They also evidenced the presence of an excess of the mid-infrared (MIR) emission in the 
LMC (70 fim excess) which is localized mostly in the atomic medium and in a large circular region around 
30-Doradus . They modeled the excess using a modified size distribution for the VSGs and interpreted it as 
a possible sign of shattering of larger dust grains into smaller VSGs. 
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Even if we can interpret infrared color variations as abundance variations, we do not know in general 
what the origin of these variations is, nor where the different types of grains are formed. Grains could be 
destroyed in supernova shocks, but grain gro wth probably takes place in dense molecular clou ds through 
grain coagulation. It is widely accepted (e.g., AUamandola et al. 1989 : Matsuura et al. lEoOsI) . that PAH 
molecules could be formed in the expanding envelopes of carbon-rich stars. The objective of our work is to 
trace dust abundance variations across the LMC in order to gather information about the dust formation and 
evolution in the ISM. 

Sections |2] and [3] briefly summarize the observations and the processing we applied to the Spitzer 
images. In Section |4] we use an empirical method to study the infrared color ratios. In Section [5l we focus 
on the dust abundance variations across the LMC using a dust emission model. Sections |6]and|7]are devoted 
to discussion and conclusion. 



2. LMC observations 



2.1. Spitzer Space Telescope Observations 



Th e Spitzer Space Tel escope (IWemer et al. carries three instruments: the Infrared Array Camera 

(IRAC. Fazio et al. 1 12004|) which pro vides imaging in th e near and MIR domain (3.6, 4.5, 5.8 and 8.0 
yum), the Infrared Spectrograph (IRS. iHouck et al. 112004) which provides both high- and lo w-resolution 



spect roscopy at MIR wavelengths, and the Multiband Imaging Photometer (MIPS) for Spitzer (IRieke et al. 



20041) providing imaging and limited spectroscopic capabilities at three FIR wavelengths (24, 70 and 160 
fim). IRAC data have an angular resolution of 1"6 - 1"9 from 3.6 to 8 fiva, and MIPS data have an angular 
resolution of 6", 18", and 40" at 24, 70, and 160 //m respectively. 

The data were obtained as par t of the Surveying the Agents of a Galaxy's Evolution (SAGE) Spitzer 
legacy survey (iMeixner et al. The survey covered a ~ 7 x 7 deg region centered on the LMC with all 

of the IRAC and MIPS bands. We used mosaics which combined both epochs of the data into a final map 
at a resolution of 3"6 pixer^ for IRAC, and 2.49, 4.8, and 15"6 pixel"' for the MIPS 24, 70, and 160 /urn 
maps, respectively. 



2.2. IRIS data 



The IRIS data are the improved version of lat est IRAS images, corresponding to the second generation 
processing (or IRAS Sky Survey Atlas, ISSA), see lMiville-Deschenes & Lagache I (120051) . These maps have 
been reprocessed to improve sensitivity and absolute calibration. Compared to the latest version of the IRAS 
data, this new product named IRIS (for Improved Reprocessing of the IRAS Survey) benefits from a better 
zodiacal light subtraction, calibration and zero level adjustments to match the DIRBE data on large angular 
scales. With a resolution near 4' the IRIS survey covers 98% of the sky at four wavelengths (12, 25, 60 and 
100 fim). 
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3. Image Processing 



We bri efly describe here the processing we applied to each map. Details of the data treatment used here 



Bernard et al. I (120081) . All data were projected on a common grid centered at Q'2000 = OS'' 18" 



are given in 

and ^2000 - -68°42'00", with a pixel size of 2' and then smoothed to the IRIS resolution of 4'. We 
subtracted the Galactic foreground emission in order to avoid any confusion betwee n the emission frorn 
the LMC and from our G alaxy. We used the SED of high latitude galactic emission of IPwek et al. I (|l997h 
(see lBernard et al. Il2008h. and the H I 21 cm emission in the heliocentric velocity range -64 < Vhei < 100 
km s"' (|Staveley-Smith et al. Il2003h . In order to eliminate residual zodiacal emission and cosmic infrared 
background (CIB), we have subtracted a background from all images corresponding to the median intensity 
in the region located 4° from the center of the maps. Point sources in archival catalogs (Two Micron All Sky 
Survey (2MASS) catalog for IRAC and SAGE catalog for MIPS) were masked from the IRAC and MIPS 
maps before degrading their resolution to 4', to match that of the IRIS maps. For the IRIS point-source 
subtraction we used the SAGE catalogs; specifically the 160, 70, 24 and 8 //m catalogs for the IRIS 100, 60, 
25 and 12 fim maps respectively. The Spitzer and IRIS images at the 4' angular resolution are presented in 
Figures [U and HI 



4. Infrared colors 



4.1. Color-Color diagrams 



We first examine the trends in the LMC color-color diagrams (see Figure O. The Ig/Iieo and I24/I160 
ratios can be considered as proxies for the PAH and VSG relative abundances with respect to that of the 
big grains (BGs) respectively. The Iioo/Ii6() ratio is directly linked to BG equilibrium temperature and 
therefore to the interstellar radiation field (ISRF) intensity. The Ig/Iieo and I24/I16O ratios show a trend with 
the I100/I160 ratio, such that the ratios increase with the BG temperature. However, in the first case the 
dispersion is quite large, whereas in the second case the trend is more pronounced. We can conclude that the 
PAH abundance does not depend strongly on the radiation field intensity. The VSG abundance, however, 
seems to be increased in warm regions of the LMC. 

To connect the obs erved brightnesses to dust abundances more precisely, we use an updated version of 



the iDesert et al. I (119901) model, which includes an improved description of the PAH IR emis sion bands at 
3.3, 6. 2, 7.7, 8.6, 11.3, 12.7 an d 16.4 lum. This version also takes into account the results of iFlagey et al. 



( 2006 ). based on the work by Rapacioli et al. I (|2005f) including the cross sections of neutral and ionized 
PAHs. We assume the same ISRF spectral shape as in the solar neighborhood. We show in Figure |4] the 
variations of the predicted I24/I16O and Ig/Iieo with I70/I16O for different values of ISRF intensity (Xisrf), 
from 0.2 to about 42, and various VSG and PAH mass abundances relative to hydrogen. The reference MW 
dust abundances from IPesert et al. I(ll990h are Ybg = 6.4 x 10"^ YysG = 4.7 x lO""*, Ypah = 4.3 x 10^"*, 
where Yx is the mass of dust species X relative to H. In Figure HI we allowed the Yvsg and Ypah parameters 
to vary from 5 xlO"^ to 0.5. Comparing the values of the Ig/Iieo^ I24/I16O and I70/I16O ratios given in Tabled] 
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to the predicted values deduced from the model, we deduce from Figure|4]that the intensity of the ISRF in 
the diffuse LMC medium is about 4.5 < Xjsrf < 6.5, taking into account the errors on the brightness ratios. 
Similarly, the values of the VSG and PAH mass abundance are in the range 3.5 10~^ < YysG < 4.7 lO"'^ and 
3.5 lO""^ < YpAH < 4.3 10"'* respectively. The diffuse ISM in the LMC is therefore subjected to an ISRF 
which is moderately stronger than that in the solar neighborhood, and appears to contain a similar amount of 
VSGs and PAHs with respect to large grains, once the effect of heating has been properly taken into account. 



4.2. Infrared correlations 



In order to measure the color ratios, we performed linear correlations between the IR brightnesses over 
the whole LMC, using 

I(Ai) = aiI(A2) + bi (1) 



7(^2) - a2l{Ai) + b2 



(2) 



where ai and a2 are the color ratio parameters to be determined and bi and bi are constants corresponding 
to eventual residual offsets in the images. To avoid any bias in the results depending on the order of the 
axis we consider, we take the mean value between ai and 1 /a2 and we define the uncertainty on each ratio 
as the difference between these two values. In order to identify outliers in the fit, we performed iterations 
of the linear regression and removed pixels distant by more than 3-cr from the linear fit, where cr is the 
standard deviation of the distance between the data and the fit. This process reached convergence in about 
15 iterations. The results of the regression are given in Table [T] A graphical representation of the table is 
given in Figured We compare the ratios obtained in the LMC with other ratios calculated in the SMC, 
in the Galaxy, and in other galaxies. The ratios from [Bernard et al^l (2008) have been deduced from the 
total SED in the overall LMC (see their Tables 5 and 7). Therefore in this case, individual HII regions such 
as the 30-Doradus region are included, which contribute a large fra ction of the to t al brig htness. For that 
reason our results presented here are different from those obtained by [Bernard et al. ( 2008 ). The correlation 
method proposed here is more appropriate to study the diffuse ISM, in the sense that pixels dominated by 
the diffuse medium are given more importance than the ones dominated by individual HII reg ions, which are 
disregarded as outliers. The quoted color ratios in the SMC are taken from lBot et al. I (120041) and have been 
deduced from Gaussian fits of the histograms of the ratios. We find the Iioo/Iieo and leo/Iieo ratios to be 66% 
and 62% lower, respectively, in the LMC, compared to the Iioo/Ino and Igo/Iivo ratios in the SMC. These 
differences must be due to the average equilibrium temperature of the large grains in the SMC being higher 
than in the LMC. However, the I25/I160 and I12/I160 ratios, which trace mostly the VSG and PAH abundance 
relative t o BGs and which to fir st order do not depend on the BG temperature, are almost the same in the two 
galaxies. iBolatto et al. I (|2006r) explored the southwest region of the bar in the SMC an d found an overall 



Is /I74 ratio close to 0.36, lower by a factor of 4 compared to our results in the LMC. Eng elbracht et al. 



(|2005h analyzed this ratio as a function of the metallicity f o r a sa mple of 34 galaxies. They show that the 
Ig /I24 ratio decreases with metallicity. iDwek & Galliano I (l2009h observed the same trend analyzing the 
emission of several nearby galaxies. They interpreted this as resulting from a lower PAH abundance in 
low-metallicity environments. That could explain why this ratio is smaller in the SMC, which metallicity is 
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smaller than that of the LMC. According to Engelbracht et al. I (l2005h . the mean ratio below one-third solar 
metallicity is 0.08 ± 0.04, using a solar metallicity reference of 12+log(0/H)=8.7 (lAUende Prieto et al. 



20011). For galaxies with metal licity comparable to that of the LMC (between 1/3 and 1/2 solar), we expect 



from (lEngelbracht et al. 1120051) a I8/I24 ratio close to 0.4. This is significantly lower than the median value 
of 1.6 we derived here for the diffuse emission in the LMC. The PAH relative abundance therefore seems 
to be larger in the di ffuse ISM of the LMC t han in other galaxies with similar low metallicity. However, 



Bernard et al. I (|2008h found a similar ratio as Engelbracht et al. I (|2005h . which could indicate that the ratio 



in external galaxies is affected by the contribution from major HIT region similar to 30-Doradus in the LMC. 



The I5.8/I8 ratio has been studied by lArendt et al. I (120081) in the Galactic center and by iFlagey et al. 



teOOd) in the Galactic diffuse ISM, as observed within the Gal actic First Look Survey ( GFLS) and Galactic 
Legacy Infrared Mid-Plan Survey E xtraordinaire (GLIMPSE. IChurchwell et al. Il2005|) . In the case of the 
Galactic center. 



Arendt et al. 



(|2008r) found an uniform ratio for the diffuse emission, independent of the 8 
jum brightness. T heir median value o f Ij.g/Is-O.S? ± 0.03 is close to our ratio of 0.33 ± 0.2. The IRAC 
ratio obtained by iFlagey et al. I (|2006r) ranges from 0.26 to 0.37 and is also in agreement with our result. 
Therefore, the h.g/l^ ratio seems to be rather similar, around 0.3, in the LMC and the MW, independent of 
metallicity. 



We find a I4.5/I8 ratio in the LMC equal to 0.18 ±0.37. iFlagey et al. I (|2006h found values between 
3.7x10"'^ and 6.5x10"^ in the MW, lower than our result by a factor of 2.8 and 4.9 respectively. The model 
we use shows that this ratio is rather independent of the ionization state of the PAHs, increasing by only 
about 10% when going from fully neutral to fully ionized PAHs. It is therefore likely that ionization is not 
the reason for the observed difference. However the size distribution of PAHs is likely to affect this ratio, 
leading to larger ratio values for smaller PAH sizes. Therefore, the observed difference could indicate the 
presence of smaller PAHs in the LMC than in the MW. 

For the I^ 6 /Is ratio we find in t he LMC a value of (2.9 ± 6.3)xlO~^, two to three times lower than the 
ratio derived by IFlagey et al. I (|2006r) (after point-source subtraction). This ratio is sensitive to PAH ioniza- 
tion and decreases by 90% when going from fully neutral to fully ionized PAHs. Ther efore, the observed 



difference with the MW value could indicate the presence of ionized PAHs in the LMC . lEngelbracht et al. 



(|2005h found that in subsolar metallicity galaxies, galaxies with high values of I4.5/I8 (> 0.1) tend to have 
low values of I8/I24 (essentially below 0.1). In the LMC we have both large I8/I24 and large I4.5/I8 ratios 
and a low metallicity compared to the solar value. The differences could be due to distinct properties of 
the near-infrared continuum or of the PAHs. Nevertheless, the emission at 3.6 yum is dominated by stellar 
emission, and therefore we caution that point-source subtraction for these data is difficult and could have an 
impact on the results. 
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5. Dust abundance variations 



5.1. Method 



We consider the IR brightness as the sum of the contributions from the 3 dust components, BG, VSG, 
and PAH : 

Iy{A) = 7^^7f ^(i) + Y^'^I^'^iA) + Y''^"C{A) (3) 

where Ybg» Yvsg and Ypah correspond to the mass abundances for each du st component . The l y values for 
each dust component in equation [3] are computed using the version of the Desert et al. I (Il990l) model de- 
scribed in Section l4Tn The intensities are computed in the photometric channel of the instruments considered 
in the study and have been corrected for the color correction factor, computed using the filter transmission 
and flux convention of each instrument used. The mass abundances are determined by fitting to the data. 



Bernard et al. I (120081) have shown that a significant emission exces s is presen t at 70 fim in the LMC 



(see their Section 5.3). This excess has also been found in t he SMC bv iBot et al. I (12004) . In the LMC, it 
is apparently associated with the neutr al phase of the ISM (Bernard et al. Il2008h and can be attributed to 
regions of enhanced VSG abundance. [Bernard et al. I (|2008l ) studied difi'erent possible origins to explain 
the 70 fim excess. The most likely explanation is variations of the VSG size distribution, and the excess 
indicates a greater abundance of large VSGs than in our Galaxy. In the context of the model we use, the 
VSG size distribution is given by 

^" oc a-"^"' (4) 



da 

where dn is the numerical dens ity of grains v yith a dimension between a and a+da, and gysn is the VSG 
size distribution. In the Galaxy. iDesert et al. I (Il990l) derived an average value of 2.6. [Bernard et al. I (120081) 
found that a value of aysG ~ 1 is required to explain the LMC observations. 



Using the 



Desert et al. 



(119901) model, we can study the variations of the different brightness ratios as a 
function of the Xjsrf and ayso- Figure[6]shows the evolution of the three ratios IieoAioo. IieoAvo and I160/I24 
for different values of the Xisrf and aysG derived using the model. The figure shows that avsG does not 
significantly affect the I160/I100 and I160/I70 ratios for Xjsrf values larger than about 0.1. For our study, we 
do not have enough spectral coverage to warrant extra free parameters, so we decided to consider only two 
values of avsG : the Galactic value (q'vsg-2.6) and the value used in Bernard et al. (2008; avsG=l)- 

As a first step we computed the normalized brightnesses between 160 /jm and 8 //m (below this wave- 
length the presence of an additional NIR continuum is generally needed to explain the obser vations) for a 



Desert et al. 



grid o f Xisrf values and the two selected values of avsG> using the MW dust abundances from 
(|l990l ) (see Section 14.11) . The brightness for each dust component i can be written ly = L^(avsG> Xisrf). 
We obtain a table of brightnesses for the R\Hs, VSGs, and BGs separately. We compare the brightness of 
our infrared maps with the values in the table at each pixel of the map, in order to derive maps of the mass 
abundances. To ensure that the Xisrf value is calculated correctly, we forced the model to reproduce the 
data at 100 and 160 /urn by increasing the weight of these data points in the SEDs. We derived the mass 
abundances minimizing the sum of the squared difference between the model and observations (y^) for each 
value of Xisrf and ayso and selected the solution with minimum 
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One shortcoming of our modeling is the assumption of a single Xbrf and therefore BG equilibrium 
temperature, along any line of sight. This is expected to be a good approximation in the LMC due to its 
near face-on viewing geometry. To account for the possibility o f a mixture of the r adiation field intensities 
along the line of sight, we follow the description proposed by bale et al. ( 2001 ) assuming a power-law 
distribution of the dust mass M(i(Xisrf) subjected to a given ISRF intensity on a given line of sight 



dMdiXisRF) ^ XJs^jfpdXisRF 0.3 < XisRF < 10^ 



(5) 



where o- is a parameter that represents the relative co ntribution of radiatio n field intensities, in the range 
a= 1-2.5. We first compute emission spectra using the iDesert et al. I (119901) model at various Xjsrf values 
(noted I™°'^(Xisrf)), then sum them over Xjsrf according to 



jtot _ 



jX,srf ^ISRF 



(6) 



We follow the same methodology as described previously, constructing a table of norm alized bright 
nesses for different a values for each grain component, using the MW abundance values from lDesert et al. 
(|l990l ) as a reference. We use the same values of aysG as derived previously and we minimize in each 
pixel with respect to the above table, determining a. The abundances of the various dust components are 
then obtained in the same way as previously. 



5.2. Results 

The objective of our modeling is to determine the distribution of the VSG and PAH abundance com- 
pared to the BG abundance across the LMC. In order to avoid any confusion between the dust emission and 
unresolved stars in the 8 /urn map which could have been left after the initial 8 /im point-source subtraction, 
we also masked the point sources from the initial 3.6 //m map, and smoothed the map to 4'. The median 
ratio between the point-source fluxes at 3.6 //m and at 8 /urn is 4. In order to remove potential contamination 
by unresolved stars at 8 jum, we scaled the 3.6 //m map by a factor 1/4 and subtracted it from the 8 //m map, 
assuming negligible ISM emission at 3.6 //m. This is an important step, specially in the stellar bar where 
the star density is high and where we identify a higher PAH abundance derived from the 8 /im emission map 
(see Section [5. 2. II) . 

Figures 17] and |9] show the spatial distribution of the PAH and VSG abundances relative to BGs. It can 
be seen that the distributions of these two dust components are quite different. The white regions around 
30-Doradus have very high;^^ values and have been masked in all the maps. In the following sections, we 
analyze separately the variations of the relative abundances of each dust component. 
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5.2.1. PAH distribution 

Figure|7] shows that the PAH relative abundance is larger than average over a large elongated region 
toward the center of the LMC. This region is located within the old-population stellar bar, delineated in 
the figure by the 2MASS contours. The PAH relative abundance increases by a factor of ~2 inside this 
region, compared to the rest of the LMC. Figure [8] shows that the PA H distribution is also s ometimes well 



correlated with the molecular clouds. This result is in agreement with 'Bernard et al. I (|2008h who proposed 
that this may be due to the presence of enriched PAH halos similar to those found in the solar neighborhood, 
mostly unresolved at the distance of the LMC. Our results also show increased PAH abundance toward some 
HII regions. This could result from the fact that HII regions are generally associated to molecular clouds. 
We also note that HII regions with high Xjsrf show low values of the PAH abundance, in agreement with 
findings in the MW that PAHs are deficient in warm HII regions. Note also the presence in Figure |7] of 
a north-south elongated region between 30-Doradus and the bar (around a = 5*^27™) with lower apparent 



PAH abundance. This is likely an artifact of the method used. This region was shown in ^ Bernard et al. 



(l2008i) to have an increased BG temperature, which they attributed to an artifact in the IRAS data at 100 fxra 



following crossing of the bright 30-Doradus region along the IRAS scans. This results in an overestimated 
ISRF and an underestimated PAH abundance. A similar effect is seen on Figure |9]for the VSG abundance. 
We also investigated whether the increase of the PAH relative abundance in the stellar bar could be due to 
a difference in the ISRF spectrum compared to what we used in the model. Indeed, the ISRF in the bar is 
likely to be dominated by the many old-population stars and should be redder than in the rest of the LMC. 
Since BGs are heated equally by UV and VIS photons, our procedure to derive the radiation field intensity 
from their temperature would lead to overestimating the UV part of the ISRF in that case. Since PAHs are 
expected to be mostly sensitive to UV photons, this in turn would lead to underestimating their abundance. 
We, however observe an opposite behavior. We therefore can exclude that the increased PAH abundance in 
the stellar bar could be due to this effect. 

The top panel of Figure [lO] shows the median SED, computed in regions 1 and 2 (see Figure ItJi with 
low and high PAH relative abundance (Ypah/Ybg = 0.076 and 0. 16), respectively. The SEDs are compared 
to the prediction of the model described by Equation [3] The error bars correspond to the standard deviation 
in the background region (same background region as described in Section O . We took into account the 
point-source subtraction in the error estimate by dividing the standard deviations by a mask whose values 
ranged between and 1 for pixels with or without contamination by point sources, respectively. Table |2] 
summarizes the coordinates, size, brightness values, and model predictions for each region. We can see that 
the 160 fim brightness for the two regions are within 15% of each other, while the difference in brightness 
at 8 jum is a factor of 2, which illustrates why the model fits give Ypah/Ybg higher by a factor of 2 for 
region 1. However, we can see that the model for the two SEDs is not able to reproduce the IRIS data at 12 
jum. The model we consider here assumes a fixed fraction of neutral and ionized PAHs of 90% and 10%, 
respectively. Therefore we suspect that the difference between the data and model at 12 fim could be due to 
this limitation of the model. We will discuss this point in Section [6l 
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5.2.2. VSG distribution 

Figure |9] shows the VSG relative abundance map. The VSG distribution across the LMC is dilferent 
from that of the PAHs. Regions with enhanced VSGs are located around 30-Doradus and in a patch near the 
center of the LMC. As for he PAHs, the increase of the VSG abundance near the center of the bar cannot 
come from a reddened ISRF, since this would produce a decrease instead of an increase of their abundance. 
As opposed to the PAH abundance, the VSG relative abundance follows the distribution of the HIT regions 
quite well (symbols in Figure |9ll. Note that a similar artifact as in the PAH abundance is evidenced in the 
north-south region between 30-Doradus and the bar (see Section [5.2. ll for details). 

The bottom panel of Figure [10] shows the median spectra for two regions (regions 3 and 4) with low and 
high VSG relative abundance (Yvsg/Ybg = 0.24 and 1.18) respectively and their corresponding models. 
The main difference between the two spectra in the FIR domain (from 60 to 160 fim) is the flatness of the 
spectrum in region 3 with high VSG relative abundance. The avsG value is equal to 1 for the model spectra, 
but the derived VSG relative abundance is increased by a factor of 2 from region 4 to 3. As discussed in 
Section [STTl we only allow two values for the VSG size distribution slope, avsG=l or the Galactic value, 
2.6. For most pixels, the value deduced from the fit is 1, but some regions require the Galactic value. This is 
the case of high Xjsrf regions such as 30-Doradus and some individual HIT regions. In these warm regions, 
the 70 yum emission is actually dominated by warm BG emission and little additional emission from VSGs 
is required. In those cases, the model with aysG - 2.6 sometimes overestimates the 70 /urn emission, which 
could indicate the presence of a steeper VSG size distribution than in the MW. In summary, results of the 
fits give QfysG = 1 in the overall LMC except for HII regions where the Galactic value (or higher value) is 
needed. 



5.2.3. ISRF mixture results 



The a pararn eter (radiation field probability distribution) used in Equation [5] ranges from 1 to 2.5. 
Dale et al. I ([200 ih describe the two extreme cases : Qf=2.5 corresponds to diffuse regions, whereas a=l 
represents photodissociation regions near young stars. In the LMC, regions with a=2.5 are more common. 
The decrease of the parameter g acr oss the LMC follows the increase of the BG equilibrium temperature 
(See Figure 7 of bemard et al. Iboosi) . 



The PAH and VSG relative abundances derived for the case where we let the ISRF strength mixture to 
vary are given in Figures \TT\ and [121 respectively. The scale in the two maps are the same as in Figures [7] 
and [21 Figure[TT] still shows the increase of the PAH abundance in the old-population stellar bar, although 
the amplitude of the variation is slightly less than that in Figure [71 The VSG relative abundance map is quite 
similar to the one shown in Figure [9l : the VSG relative abundance still follows roughly the distribution of 
the HII regions and the main region of enhanced VSG abundance is still located around 30-Doradus and 
near the center of the LMC. Results do not appear to change significantly when we consider the effect of a 
reahstic ISRF strength mixture along the line of sight. 
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6. Discussion 



PAHs and VSGs could have a different origin or could be subjected to different processing in the ISM. 
Several studies of the MIR aromatic bands and cont inuum emission arou nd 15 fim showed that PAH and 
VSG emission are linked to star formation activity. Forster et al. I (120041) describe the PAH component as 
representative of the quiescent environments while the VSG component traces the active star formation sites. 
They attribute the disappearance of PAH bands in HII regions to the destruction of PAHs in regions with high 
intensity of the radi ation field (in agreemen t with the decrease of t he IRAS 12/100 ratio in regions of intense 
UV radiation, from lsoulanger et al. lllQSsL Povich et al. ( 2007 ) considered two possible mechanisms for 
PAH destruction in HII regions : stellar UV photons and X rays . They conclud ed that only UV photons 
could be responsible for the absence of PAHs in these regions. IContursi et all hood) examined the N4 
region in the LMC using ISOCAM data and found that the maximum emission is located at the interface 
be tween the HII region an d the molecula r cloud , as was also found from observations of M17 in our Galaxy 
by lCesarsky et al. I (119961) . iHelou et al. 1 J2OO4') studied the spiral galaxy NGC 300 using Spitzer data and 
found that emission at 24 /im traces the star-forming regions whereas the emission at 8 /urn highlights 
the edges of these regions. However, interpreting those observations in terms of VSG abundance remains 
difficult, since the excitation conditions are generally poorly constrained. 

At least a fraction of the PAH population is thought to be formed and inj ected into the ISM by the 
stellar winds following mass loss from old carbon-rich stars. iLoup et al. I (119991) conducted a survey of 0.5 
square deg^ in the bar of the LMC using the ISOCAM and DENIS instruments. Their objective was to 
build a complete sample of thermally pulsing asymptotic giant branch (TP-AGB) stars. During this phase 
stars lose mass and carbon-rich AGB could contribute significantly to the formation of PAHs. They found 
approximately 300 TP-AGB stars, of which 9% have high mass-loss rates. Some studies of post-AGB stars 
in tra nsition between the AGB and planet ar y nebulae phase e xhibit signatures of dust attributed to PAHs 
(e.g., Venkata Raman & Anandarao 2008 ). Justtanont et al. I ( 1996 ) analyzed MIR spectroscopy in three 
post-AGB stars, known to exhibit an emission feature at 21 fim. In addition to a weak emission at 1 1.3 jum, 
they found several features which match the emission bands in the spectrum of chrysene (C8H12), which is 
one of the simplest PAH molecules. The 21 /jm emission feature has never been observed in carbon rich 
AGB stars, nor planetary nebulae. The PAH molecule responsible for this feature could be f ormed in the dust 



shells around some post-AGB stars and lead to increase the abundance of this PAH species (IJusttanont et al. 
19961) . 



In this study, we showed that the PAH relative abundance in the LMC is maximum in an extended region 
corresponding to the old-population stellar bar, a region with no strong current stellar formation activity. To 
our knowledge, this is the first time such an effect is evidenced at Galactic scale. We interpret this finding as 
an indication that PAHs may have formed efficiently in this region in the past from carbon-rich stars winds 
and have remained overabundant in the surrounding ISM. Our result also shows the presence of a region 
with a significantly increased VSG abundance in the same region, although with a lesser spatial extent. This 
is hkely to reveal that VSGs also formed efficiently from the same process. The difference in the spatial 
extent of the regions could either reflect specific formation processes or differences between the evolution 
of the two species in the ISM. In the later hypothesis, our result would suggest that PAHs are in fact more 
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resistant than VSGs to destruction processes in the diffuse ISM. 



A careful examination of Figure [13] which compares the spatial distribution of PAHs and VSGs in- 
dicates that there aie other regions where the PAH and VSG abundances are spatially correlated. Most 
noticeably, the region surrounding 30-Doradus and its extension to the southwest both have enhanced PAH 
and VSG abundances. However, significant differences are also apparent. For instance, the region of en- 
hanced PAH abundance to the west of the main molecular ridge. South of 30-Doradus , as well as the arc-like 
molecular ridge to the east have no counterpart in the VSG abundance map. Similarly, at the northwest end 
of the stellar bar, the PAH and VSG abundance maps appear anticorrelated. 

In the solar neighb orhood, such independ ent variations of the distribution of PAHs and VSGs have been 
observed. For instance, Bernard et al. I (Il993h showed the existence of increased PAH abundance regions at 
the periphery of some translucent clouds. This could result from the destruction of larger species such as 
VSGs into PAH molecule s , indu ced by the action of external low-velocity shocks at the surface of the clouds. 
MiviUe-Deschenes et al. I (120021) evidenced large variations of the PAH emission in the Ursa Major cirrus. 
They proposed that the observed color variations are due to actual PAH abundance variations due to grain 
shattering in energetic grain-grain collisions induced by turbulence or grain coagulation caused by small 
turbulent motions. Therefore regions with large turbulent motion could generate PAHs through shattering 
of larger dust particles, such as VSGs. Indeed, the destru ction of graphitic d ust is expected to lead to the 



formation of large two dimensional PAH-like molecules. iBeme et al. I (|2007n studied a photodissociation 
region using the Spitzer IRS da t a and found a disappearance of VSG emission while the PAH emission 
increases. Like ICesarsky et al. I (|2000l ) they suggested the possibility that VSGs could evolve into PAHs 
under UV photon driven processes. Under this hypothesis, VSGs can be viewed as PAH clusters. Although 
at the distance of the LMC, such small-scale structure would not be resolved, it is likely that the differences 
we evidence between the spatial distribution of PAHs and VSGs are of a similar nature. It is worth noticing 
that, in areas with no VSG enhancement such as the southern tip of the arc-shaped molecular ridge, the 
PAH abundance enhancement appears to trace the surface of the molecular clouds, a similar situation as that 
observed in nearby halo clouds, although with much larger physical scale, given the distance of the LMC. 



Similar variations were observed in the SMC by iBolatto et al. I (|2006r) who found large variations of 
the Is /I24 ratio across the southwest region of the bar. They attributed these to spatial variations of the PAH 
abundance, with the PAH abundance decreasing with increasing 24 fim brightness. They suggested that the 
PAH abundance varia tions are essentially driven by the local ISRF and PAH photodestruction. In the LMC, 
Bernard et al. I (120081) showed that a large region around of 30-Doradus has excess 70 jum emission, which 
they interpreted through the flattening of the VSG size distribution, potentially caused by the erosion of 
larger BGs by shocks propagating from the 30-Doradus region. In our study, we find that this same region 
also has an increased VSG abundance. Our result shows that this region, in addition to showing a n increase 
of the largest VSGs also shows an increase of the overall VSG abundance. [Bernard et al. I (|2008|) indicated 
that the increase of the largest VSG abundance would correspond to - 13% of the BG mass, which converts 
into a VSG abundance increase of 77% for a standard MW dust mixture. The excess we measure in this 
region corresponds to an increased abundance by a factor of about 2.4, and is th erefore larger than that 
imputable to VSG size distribution change only. Similarly, iLisenfeld et al. I (120021) analyzed the emission 
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from NGC 1569, a dwarf galaxy with a metallicity between those of the LMC and the SMC. Using the 



Desert et al. I (119901) model, they evidenced an enhancement of the VSG abundance relative to BGs, with 
respect to the solar neighborhood, by a factor of 2-7, depending on the model parameters they used. They 
explain this result by the BG destruction due to shocks in the turbulent ISM of this galaxy. 

We therefore consider it likely that the differences evidenced here between the spatial distribution of the 
PAH and VSG abundances relative to BGs are tracing both the production of PAHs from the destruction of 
VSGs at the surface of molecular clouds, and the production of VSGs from the destruction of larger grains 
in shocks. 

We saw in Figure [TOl that the model with neutral PAHs was not able to reproduce the IRIS data at 12 
fim. We suspect that this is due to the ionization state of the PAHs or the nature of the molecules, such as 
dehydrogenation. It is difficult to constrain the nature of the PAHs in the LMC from the data presented here 
only. However, varying the PAH ionization fraction in the model, we find that the model requires 100% of 
the PAHs in the ionized state to explain the observations at 12 fim (see Figure [141). Even if this model does 
not reproduce the data at 12 ^m perfect l y, the difference between the model and the data is smaller than 



with our initial hypothesis. iFlagey et al. I (120061) obtained a neutral over ionized PAHs ratio close to 50 % in 
the MW (GLIMPSE field). Therefore, PAHs in the LMC could be more ionized than in the Galactic plane. 
However, the exact nature of PAHs in the LMC could be different than in our Galaxy. A study of the nature 
of the PAH molecules is essential to determine the origin of the observed difference at 12 yum between our 
Galaxy and the LMC. Such study may become possible using the data from the Spitzer Sage-Spec legacy 
program. 



7. Conclusion 

We used the data of the SAGE Spitzer legacy survey, combined to the IRIS data to measure the varia- 
tions in the colors of infrared emission from dust in the LMC. Comparisons with results from the literature 
showed a higher I8/I24 ratio in the LMC than in other low-metallicity environments, which could result from 
a relatively higher PAH abundance in the LMC with respect to its metallicity class. Morever, the I3.6/I8 and 
I4.5/I8 ratios are smaller in the LMC than in the MW. This could result from different properties of the 
near-infrared continuum or from different properties of the PAHs. In particular, PAHs could be smaller in 
the LMC, compared to our Galaxy, and more ionized. The SED modeling of regions with enhanced PAH 
and VSG relative abundances in the LMC show a similar result : the model used requires a large fraction 
(potentially 100%) of ionized PAHs in order to minimize the difference between the model and the data at 
12 /urn. 

In order to understand the processes of grain formation and destruction at the scale of the LMC, we con- 
structed for the first time maps of the PAH and VSG relative abundances with respect to BGs of the overall 
galaxy, taking into account the distribution of the radiation field intensity, as inferred from the equilibrium 
temperature of the BGs. We evidenced an increase of the PAH relative abundance over an extended region 
which matches the extent of the old-population stellar bar. PAHs are also overabundant at the periphery of 
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some molecular clouds. The VSG abundance is increased over a region around 30-Doradus and in a small 

region near the center of the bar. This clear distinction between the PAH and VSG spatial distribution could 
be due to different origins or different processings in the ISM for these two dust components. In agreement 
with several studies, the VSG abundance distribution appears to follow the star formation activity while that 
of PAHs could follow more quiescent environments. 

We interpret the larger PAH abundance in the bar as a sign of the past formation of these species by 
C-rich mass losing stars. VSGs could also be formed from the same process but the difference in the spatial 
distributions remains unclear : the two dust components could originate from different stellar populations 
with different distribution, or VSGs and PAHs could have been evolved differently in the diffuse ISM. If 
destruction is at work, this may indicate that VSGs are actually destroyed faster than PAHs in the diffuse 
ISM. We interpret the rest of the abundance variations as signs of the evolution of dust species of different 
sizes in the ISM. The region around 30-Doradus shows increased VSG abundance, which could result from 
the destruction of larger BGs in shocks. Similarly, regions close to the surface of molecular clouds with an 
increase PAH relative abundance could result from the destruction of VSGs into PAHs by shocks impacting 
the cloud's surfaces. 
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Fig. 2. — IRIS images at 100, 60, 25, and 12 /urn, at the 4' angular resolution, after the processing described 
in Section[3l The grayscale is logio(Iy(MJy sr"^)). 
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Fig. 3. — Ig/Iieo (left) and I24/I160 /^m (right) vs. I100/I160 /^m ratios. Every group of three points has 
been averaged to produce one plotted point. 




Fig. 4. — Predicted I24/I16O (left) and Ig/Iieo (right) vs. I70/I160 using the iDesert et al. I (119901) model, for 
various intensities of the ISRF (from 0.2 to 42) and VSG mass abundances (left) and PAH mass abundances 
(right) (from Yvs g pah = 5 X 10~^ to Yvsg pah = 0.5). The bold lines correspond to the reference MW dust 
abundances from lDesert et al. I (119901) : YysG = 4.7 x 10""^ and Ypah - 4.3 x 10""*. The gray rectangular 



boxes show the color values derived for the LMC (see the text). The size of the boxes shows the error. 



-20- 




Fig. 5. — SEDs deduced from the IR r atios given in Table [T] for the LMC. Our results are represented by 
circles. The SED of the whole LMC by iBemard et al. I (|2008h is shown by stars. The SMC points are from 
Bot et al. (2004; triangles) and Bolatto06 et al. (2007; cross). The Galactic diffuse ISM points are from 
Flagey06 et al. (2006; diamonds). A slight shift in wavelengths has been applied for clarity. 




Fig. 6. — I160/I100 (solid lines), I160/I70 (dot lines) and I160/I24 (dash-dotted lines) deduced from the 
Desert et al. model as a function of the intensity of the ISRF (Xjsrf) and of the index of the VSG 

size distribution (avso)- 



-21 - 




Fig. 7. — Map of the PAH relative abundance with respect to BGs at a resolution of 4', overlaid with contours 
at 0.5 and 0.7 MJy sr~^ of the Two Microns All Sky Survey (2MASS) J-band emission, convolved to the 
same resolution. The contours show the extent of the stellar bar. Regions 1 and 2 indicated by rectangles are 
discussed in Section 15.2.1 1 The horizontal line in the color bar shows the value of the standard MW relative 
abundance. 
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Fig. 8. — Map of the PAH relative abundance with respect to BGs, overlaid with the NANTEN ^^CO (J=l-0) 
integrated intensity contours convolved to the 4' resolution, at 0.4, 2 and 6 K km s~^ The hnes indicate 
the limits of the survey. The horizontal line in the color bar shows the value of the standard MW relative 
abundance. 
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Fig. 9. — Map of the VSG relative abundance with respect to BGs. The overlaid symbols show the HII 
regions (iDavies et al. 1 19761. catalog). Regions 3 and 4 indicated by rectangles are described in Table|2l The 
horizontal line in the color bar shows the value of the standard MW relative abundance. 
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Fig. 10. — Median brightnesses between 5.8 and 160 /im deduced from the data (diamond linked by contin- 
uous line) and brightnesses deduced from the model (stars linked by dashed line). The top panel corresponds 
to SEDs in regions 1 and 2 (see Figure |7]) with a low and high PAH relative abundance in black and gray 
respectively. The bottom panel corresponds to values in regions 3 and 4 (see Figure ^ with low and high 
VSG relative abundance in black and gray respectively. 
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Fig. 11. — Map of the PAH relative abundance with respect to BGs, for the case of an ISRF mixture along 
the line of sight. Contours correspond to the Two Microns All Sky Survey (2MASS) J-band map, convolved 
to the 4' resolution. The horizontal line in the color bar shows the value of the standard MW relative 
abundance. 



-26- 




Fig. 12. — Map of the VSG relative abundance with respect to BGs, for the case of an ISRF mixture along 
the line of sight. The overlaid symbols show the HII regions (Davies et al, 1976 catalog). The horizontal 
Une in the color bar shows the value of the standard MW relative abundance. 
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Fig. 13. — Map of the PAH relative abundance with respect to BGs, overlaid with contours at 0.6 and 0.9 of 
the VSG relative abundance. 
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Fig. 14. — Same figure as Figure [TOl using a model with 100% ionized PAHs. 
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Table 1 : Infrared color ratios of data. 



IR ratios 










34 galaxies'! 


T /T 

il00/il60f 


(0.0 ± V.j) IV 


A _i_ 1 n\ 1 A~l 
(/.O ± 1.9) 10 


1 A _i_ A 
l.U ± U.Z 






T^„ /T, 

^70/^160 




y^.j ± U.y ^ ru 








I60/Il60f 


(3.1 ±0.7)10-1 


(3.4 ± 0.7) 10-1 


0.5 ±0.1 






I25/Il60f 


(2.4 ± 1.0)10-2 


(1.9 ±0.6) 10-2 


(2.0 ± 1.0)10-2 






I24/I16O 


(1.7 ±0.3) 10-2 


(2.5 ± 0.5) 10-2 








Il2/Il60f 


(1.5 ±0.4) 10-2 


(7.4 ± 3.8) 10-3 


(1.1 ±0.6)10-2 






Is /1 160 


(2.5 ± 0.2) 10-2 


(1.6 ±0.4) 10-2 








I5.8/I16O 


(9.2 ±7.1) 10-3 


(2.7 ± 0.9) 10-3 








Is /I24 


1.6 ±0.3 


(6.3 + 1.4)10-1 


0.368 




(8.0 ±4.0) 10-2h 


Is/Il2 


1.7 ±0.2 


2.1 + 1.1 








Is.s/Is 


(3.3 + 2.1)10-1 


(1.7 + 0.6)10-1 




(3.7±0.3)10-i'= 
(2.6-3.7) 10-1 




I4.5/I5.8 


(5.4 ± 7.7) 10-1 


(4.3 ± 28.5) 10-2 








I3.6/I4.5 


(1.6 ±0.2) 10-1 


1.1 ±7.6 








I4.5/I8 


(1.8 ±3.7) 10-1 


(7.4 ± 5.6) 10-3 




(3.7 - 6.5) 10-2*= 


>0.1 if I8/I24 < O.l'' 


I3.6/I8 


(2.9 ± 6.3) 10-3 


(7.9 ± 7.6) 10-3 




(5.9 - 9.4) 10-3*= 





iBemard et al. 



0m (overall LMC) 



MBot et al. I (I2OO4I ) (ov erall SMC) 



Galactic Center from 



Engelbracht et al. I (l2005h 



Arendtetal. I(l2008h 



Gala ctic First Look Survey / Galactic Legacy Infrared Mid-Plane Survey Extraordinaire from lFlagey et al. 
(I2003) 



^ 170 um instead of 1 60 /urn for the lBot et al. I (120041) ratios 



Bolatto et al. 



(|2006l ) (Southwest Bar region of the SMC) 
I' In subsolar metallici ty environments (with a solar metallicity reference 
Allende Prieto et al. IboOll) 



12+log(0/H)-8.7, 
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Table 2: SEDs and Results of the SED Fits for the Four Selected Regions : 1, High 2, Low 3, 
High ^; 4, Low ^. 



YpAH . 



Region 




1 


2 


3 


4 


ffiooo 




05:23:00 


05:43:00 


05:35:00 


05:42:00 


(^2000 




-69:30:00 


-71:10:00 


-69:47:00 


-70:15:00 


Box size (°) 




0.4 X 0.4 


0.4 X 0.4 


0.4 X 0.4 


0.4 X 0.4 


Ii60 (MJy sr" 




46.72 ± 0.49 


40.30 ± 0.47 


66.05 ±0.51 


56.78 ± 0.48 


lioo (MJy sr- 




34.06 ± 0.23 


17.68 ± 0.22 


53.34 ± 0.24 


30.51 ±0.22 


I70 (MJy sr" 


') 


21.91 ±0.10 


7.12 ±0.09 


41.90 ±0.10 


18.32 ± 0.09 


leo (MJy sr" 


') 


14.67 ± 0.05 


4.23 ± 0.05 


32.70 ± 0.05 


9.38 ± 0.05 


I25 (MJy sr- 


') 


1.08 ±0.04 


(4.23 ± 0.36) 10-1 


3.20 ± 0.04 


(9.28 ± 0.36) 10 


I24 (MJy sr" 


') 


(8.02 ± 0.37) 10-1 


(3.48 ± 0.37) 10-1 


1.70 ±0.04 


(6.69 ± 0.37) 10 


I12 (MJy sr" 


') 


(9.01 ± 0.44) 10-1 


(5.22 ± 0.42) 10-1 


(9.81 ± 0.43) 10-1 


(7.34 ± 0.42) 10 


Ig (MJy sr-i 


) 


1.54 ±0.02 


(8.31 ±0.22) 10-1 


1.51 ±0.02 


1.26 ± 0.02 


Ypah/Ybg 




1.41 10-1 


7.27 10-2 


1.37 10-1 


8.31 10-2 


Yvsg/Ybg 




4.87 10-1 


1.22 10-1 


1.09 


2.28 10-1 


XlSRF 




0.85 


0.49 


0.28 


0.59 



